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Abstract. In the framework of the relativistic and gauge invariant spectral integral 
technique, we calculate radiative decays p(770) — > 77r(140) and cj(780) — > 7 7r(140) 
supposing all mesons (71", p and u) to be quark-antiquark states. The qq wave 
functions found for mesons and photon lead to a reasonably good description of data 
(r (e ± p) ± = 68 ± 30 keV, T {e xp) „ = 77 ± 28 keV, T (exp) „ = 776 ± 45 keV) that 
makes it possible to estimate the coupling for the bremsstrahlung emission of pion 
by quarks g n = gn(u — > dx). We have found two values for the pion bremsstrahlung 
coupling: \g„\ = 16.7 ± 0.3 t°i (Solution I) and \g„\ = 3.0 ± 0.3 ±%{ (Solution II). 
Within SU(6)-symmetry for nucleons, Solution I gives us for irNN coupling the value 
16.4 < gl NN /(4n) < 23.2 that is in qualitative agreement with the ttN scattering 
data, 5^jvat/(47t) ~ 14. For excited states, we have estimated the partial widths in 
Solution I as follows: T{pf s -> 7 7r) ~ 10 - 130 keV, T(p° s 7 tt) ~ 10 - 130 keV, 
T(uj2S 7 7r ) — 60 — 1080 keV. The large uncertainties emphasise the necessity to 
carry out measurements of the meson radiative processes in the region of large masses. 
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1. Introduction 

The radiative decay amplitude is a necessary element for the study of the quark-gluon 
structure of hadrons. In this paper, we present the calculation of the radiative decays of 
quark-ant iquark states (qq)i n — P, w into 771-. In this way, we continue the calculations 
initiated in pQ where radiative transitions of quarkonium states (QQ)i n — > j(QQ)out 
were studied, with the production of massive outgoing states (QQ)out- Considering 
the production of the 77T system, a particular necessity is to take into account, together 
with the annihilation qq — > tt, an additional process of the bremsstrahlung type, namely, 
q qir. 

We treat the meson decay amplitude as triangle diagram of constituent quarks 
(additive quark model) calculated in terms of the spectral integration technique, see 
[2] and references therein. The spectral integral technique is rather profitable for the 
description of composite particles, for the content of a composite system is thus strictly 
controlled. Besides, this technique is rather convenient for the description of high spin 
states. 

The equation for the composite qq systems in the spectral integration technique 
was suggested in [3], it is a direct generalisation of the dispersion N/D equation [I] 
when the iV-function was represented as an infinite sum of separable vertices, see [2] for 
detail. In terms of this equation, the bb and cc quarkonia were considered in [5], while 
the light-quark qq mesons were studied in [6] . 

In [6], the levels of the one-component qq systems (with / = 1 or / = which are 
almost pure ss or nn = (uu + dd) / \/2 states) were reconstructed as well as their wave 
functions. The qq systems are formed at distances, where perturbative QCD does not 
work (r ~ 0.5 — 1.0 fm). In this region (the region of soft interactions), we deal with 
constituent quarks and effective massive gluons (with mass of the order of 700-1000 MeV 
[3 El El QUI EI])- It means that quark-ant iquark interactions undergone a significant 
changes as compared to small distances; besides, at large distances the confinement 
forces work. Therefore, interactions in the soft region should be reconstructed on the 
basis of experimental studies - in [6], the qq interaction was reconstructed on the basis 
of available data for gg-levels and the gg-meson radiative decays. 

The standard way to investigate quark-antiquark systems is to apply the Bethe- 
Salpeter equation [12] written in terms of Feynman integrals. One may find the examples 
of such a study of light quark-antiquark systems in [131 UU EE EE Ej and for heavy 
quarkonia (cc and bb) in [T71 [HI [191 12Q1 [211 122], see also references therein. 

However, one should keep in mind an important difference between the standard 
Bethe-Salpeter equation and that written in terms of the spectral integral [3j. In the 
dispersion relation technique, the constituents in the intermediate state are mass-on- 
shell, kf = m 2 , while in the Feynman technique, which is used in the Bethe-Salpeter 
equation, kf 7^ m 2 . So, in the spectral integral equation, when the high spin state 
structures are calculated, we have a numerical factor kf = m 2 , while in the Feynman 
technique one should write kf = m 2 + (kf — m 2 ). Here, the first term in the right- 
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hand side provides us the contribution similar to that used in the spectral integration 
technique, while the second term cancels one of denominators of the kernel of the Bethe- 
Salpeter equation, that results in the penguin or tadpole type diagrams - let us call them 
zoo-diagrams. A particular property of the spectral integral technique is the exclusion 
of zoo-diagrams from the equation for composite systems. 

The spectral integral equation [3] gives us a unique solution for the quark-ant iquark 
levels and their wave functions, provided the interquark interaction is known. Let 
us emphasize that the equation works for both instantaneous interactions and the t- 
channel exchanges with retardation, and even for the energy-dependent interactions: this 
follows from the fact that the equation itself is the modified dispersion relation for the 
amplitude. For solving the inverse problem, that is, for reconstructing the interaction, 
it is not enough to know the meson masses — one should know wave functions of quark- 
antiquark systems. Such an information is contained in the hadronic form factors and 
radiative decay amplitudes. Therefore, in the approach of refs. [21 El [6] , we consider 
simultaneously the meson levels in terms of the spectral integral equations and the 
meson radiative transitions in terms of the double dispersion relations over qq states (or 
over corresponding meson masses) — in this way all calculations are carried out within 
compatible methods. 

The calculation of radiative transition amplitudes in terms of the double dispersive 
integrals was performed for some selected reactions in [231 EU ESI [261 [27] - the basic 
points of the method of operator expansion used in the calculation of double dispersive 
integrals can be found in [U [2j [28] . 

The analyses of the light qq systems [6] and heavy QQ quarkonia [5] in terms of 
the spectral integral equation differ from one another in certain respect, because the 
available experimental data are of different sort: for the QQ systems the only known 
are low-lying states (with an exception for the 1 quarkonia T and ip where a long 
series of vector states was discovered in the e + e~ annihilation). At the same time, for the 
low-lying states there exists a rich set of data on radiative decays: (QQ)m — > l(QQ)out 
and (QQ) m — > 77. For the light quark sector (qq systems), there exists an abundant 
information on masses of highly excited states with different J PC (see [291 ESS EE E21 [33] 
and surveys (2j [3H [35]), but the knowledge of radiative decays is rather poor. 

Despite the scarcity of data on radiative decays, the light qq states have been studied 
in [6], relying upon our knowledge of linear trajectories in the (n, M 2 )-plane, where n is 
the radial quantum number of the qq- meson with mass M (see [2, [36]). We hope that it 
may somehow compensate the lack of information on the wave functions. In the fitting 
procedure [6] , the main attention was paid to the states with large masses, expecting to 
extract the confinement interaction. We obtained that the strong t-channel interaction 
(which, as we think, determines the confinement) should exist in both scalar / <g> I and 
vector 7^ ® 7^ channels. The fitting results point rather reliably to the equality of these 
t-channel interactions [6J. 

Obviously, the fitting results presented in [6] should be checked (and, if necessary, 
improved) by investigating the other radiative decays - following to this program we 
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consider here the decays p — > 77? and 00 — > 77?. Small mass of the pion requires to take 
into account not only the process of photon emission with a subsequent quark-ant iquark 
annihilation qq — > n (triangle diagram of the additive quark model, Fig. CO) but also 
the bremsstrahlung-type emission of pion q — > irq, with subsequent quark-ant iquark 
annihilation into photon qq — > 7, see Fig. El Therefore, the key points in the calculation 
of the p, to — > 77T decays is to know qq wave functions of pion and vector mesons (p and 
uj) as well as qq wave function of the photon 7 — > qq. Also the fitting procedure calls us 
to determine the pion bremsstrahlung constant for the process q — > nq. 




Figure 1. a), c) Triangle diagrams for radiative transition (qq)i„ — > 77r with the 
emission of photon by quark; here p 2 = Mf n , p' 2 = M 2 and (p — p') 2 = q 2 . b), d) 
Cuttings of the triangle diagrams []Jl and [Tfc signify the double discontinuity of the 
spectral integral with intermediate-state momentum squared P 2 = s, P' 2 = s' and 
{P-P') 2 = q 2 . 

1.1. Photon wave function 

For the region < Q 2 < 1 (GeV/c) 2 (here Q 2 = — q 2 ), the light-quark components of 
the photon wave function 7*(Q 2 ) — > 99 (<? = u,d,s) are determined in [37] (see also 
[2]) on the basis of data for the transitions 7[ ,i],i]' — > 77*(Q 2 ) and reactions of e + e _ - 
annihilation: e + e~ — * p°,u,(j) and e + e~ — > hadrons at 1 < E e + e - < 3.7 GeV (in a more 
rough approximation the wave function 7(Q 2 ) — > qq was found in [38]). 

Conventionally, one may consider two pieces of the photon wave function: soft and 
hard ones. Hard component relates to the point-like vertex 7 — > qq, it is responsible 
for the production of quark-antiquark pair at high virtuality. At high energies of the 
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Figure 2. a), c) Triangle diagrams for radiative transition (qq)i n — > 77r with the 
emission of pion by quark; here p 2 = M 2 n , p 2 = M%, b), d) Cuttings of the triangle 
diagramsEk and[2fc for getting double discontinuity of the spectral integral with P 2 = s, 
P' 2 = s' and {P-P'f =p\. 

e + e~ system, the ratio of cross sections R = o~(e + e~ — > hadrons) / a(e + e~ — > is 
determined by the hard component of photon wave function, while soft component is 
responsible for the production of low-energy quark-antiquark vector states such as p°, 
u, 0(1020), and their excitations. 

In the spectral integral technique, the quark wave function of the photon, 7*(Q 2 ) — ► 
qq, is defined as follows: 

MoW*) = . 

where G^ q g(s) is the vertex for the transition of photon into qq state, depending on 
the invariant energy squared, s, of qq system. In terms of the light-cone variables 
s = (m 2 + k±)/[x(l — x)], where m is the quark mass, k± and x are the light-cone 
characteristics of quarks: transverse momentum and a part of longitudinal momentum. 

Rather schematically, the vertex function G? 7 _ >g g(s) may be divided into two terms. 
The first term is responsible for the soft component which is due to the transition of 
photon to vector qq meson 7 —>■ V —>■ qq, while the second one describes the point-like 
interaction in the hard domain. The principal characteristics of the soft component is 
the threshold value of the vertex and the rate of its decrease with energy. The hard 
component of the vertex is characterized by the energy where the point-like interaction 
becomes dominant. 
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In [38j, the photon wave function has been found assuming the quark relative 
momentum dependence to be the same for all quark vertices: g^uuik 2 ) = g^ddi^ 2 ) 
= g-y-* S s{k 2 ), where we redenoted Gj^ q q(s) — > g 1 -^ q q(k 2 ) with k 2 = s/A — m 2 . The 
hypothesis of the vertex universality for u and d quarks used in [37] . 

Gj^ U u(s) = G 7 _> d( j(s) = G 7 (s) , (2) 

looks rather trustworthy because of the degeneracy of p and uj states, though the 
similarity in the fc-dependence for non-strange and strange quarks may be violated. 
Using experimental data on the transitions 77*(Q 2 ) — > 7r°, ?7, 77' only, one cannot 
determine the parameters (C,b,s - see below Eqs. fl3J) and PI) for both G 7 ^ sf (s) 
and G 7 (s). We also add the e + e~ annihilation data for the determination of wave 
functions, that is e + e~ — > 7* — > p°, uj, 0(1020), together with the ratio R(E e + e -) = 
cr(e + e~ — > hadrons) / a(e + e~ — > p + n~) at E e + e - > 1 GeV. The reactions e + e~ — > 7* — > 
p°, uj, 0(1020) are rather sensitive to the parameters C a ,b a , while the data on R(E e + e -) 
allow us to fix the parameter s - 

The transition vertices for uu, dd — > 7 have been chosen in the form: 

uu, dd: G 7 (s) = c 7 (^exp(-bjs) + c\ exp(-&]s)^ + - - e \ {s _ s ^ ( 3 ) 

and the following parameter values have been found [21 [37j : 

mm, dd : c 7 = 32.506, c 7 , = -0.0187, bj = 4 GeV" 2 , 6] = 0.8GeV~ 2 , (4) 
62 = 15 GeV" 2 , s 7 , = 1.614 GeV 2 . 

With these parameters, we have a good description of the available experimental data 
for V — > e + e~ and two-photon decays, see [21 E]. 

TTie p, uj and tt wave functions 
We characterise gg-states by the following momentum-dependent wave functions: 

^ j \k 2 ) = c rSx2 > (5) 



where S, L, J are the spin, orbital momentum and total momentum of the qq system 
with mass Mi s ' L ' J \ 



1.2.1. p{nL) andujinV) states We introduce spin-orbital operators and wave functions 
for the states with dominant L = 0, 2 as follows: 



L = 
dominant L 







dominant L = 2 



r 



^(0,0,0)^2) 

7^ W) (fc 2 ) 



3/V2 



(6) 



^ 1A1) (* 2 ) 



Here k~ 



is the relative quark-antiquark momentum, kj^ 
so fc- 1 _L P = ki + k 2 ; likewise, ~ 1 



(g^ - P^/P 2 )^ = 

iii' = so k 1 - 1 P = ki + k 2 ; likewise, 7^ = 5^7^. Definition of spin- 

momentum operators for other states can be found in [21 [2H] ■ 
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Generally, the states with different L mix with each other: 

^(M) (s) = C£V^ 1A1) (* 2 ) + (k^ - \khi) ^\k%{7) 



But, according to [0], we have with a good accuracy = = 0, so below we put 



^10 — °22 — L - 



We parameterise the qq wave functions of p, uj states, ^^^'^{k 2 ), with the following 



formula: 

li 



<) L,J) (^ 2 ) =e-^ 2 Y,c t (S,L,J ] n)\kr 1 , (8) 



i=l 



with cutting parameter /3 = 1.2 GeV~ 2 . In Eq. ([8]), we use the notation \k\ — \Js/4 — m 2 
(m is the mass of the light constituent quark, m ~ 350 MeV). 

The constants Ci(S, L, J;n), in GeV units, for mesons with L = 0, ^ ,s ' i=0 ' J ^(/c 2 ), 
and L = 2, ^ 5 ' L=2 ' J ^(fc 2 ), are presented in Eq. (9) and (10). 



2 




Figure 3. Trajectory for p n s and u n s states found in [5] (M p („s) = A/ W („s)). 
Experimental values of the masses on p- and w-trajectories are equal to: [pi5(775± 10), 
p 2S (1460 ± 20), p 3 s(1870 ± 70), /j 4 s(2H0 ± 35)] and [w is (782), w 2S (1430 ± 50), 
w 3 s(~ 1830), w 4 s(2205 ± 40)]. 



In the solution found in [6], the p„5 and uj n s mesons are degenerated: M p ( n s) = 
Mu(nS)j see Fig- El Coefficients Ci(S — 1, L — 0, J = 1; n) for n < 4 (recall that n is 
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radial excitation number) read: 





p(lS),u(lS) 


p(2S),u(2S) 


p(3S),u(3S) 


p(AS),u(AS) 


i 




/1 nil 

^" j 
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44.2 


-47.0 


34.4 


256.1 
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147.9 


96.4 


367.3 


-3816.4 
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-2576.7 


1694.4 


-6627.1 


21285.8 


4 


10145.9 


-8835.1 


31300.6 


-61891.6 


5 


-20331.5 


18954.3 


-72495.7 


106967.9 


6 


23805.7 


-21715.0 


95497.7 


-115547.6 


7 


-16569.8 


13585.9 


-73882.6 


77608.2 


8 


6338.4 


-3952.2 


31633.5 


-29980.2 


9 


-941.1 


119.3 


-5588.5 


4927.5 


10 


-59.0 


26.4 


-333.1 


258.1 


11 


-16.0 


88.7 


43.2 


-25.9 



(9) 



For p and mesons with dominant L = 2 we have the following Cj(5' = 1,L = 2, J = 
l;n): 





P (lD),co{lD) 


p{2D),u{2D) 


p{3D),u{3D) 


p(4D),u(4D) 




1 (1,2,1) 




1 (1-2,1) 




1 


32.6 


1.9 


295.8 


1109.3 


2 


-297.9 


-20.8 


-2587.2 


-9686.9 


3 


1030.3 


85.0 


8635.8 


32404.0 


4 


-1720.3 


-207.3 


-13721.7 


-52043.5 


5 


1257.2 


242.8 


9530.7 


36934.5 


6 


68.1 


4.0 


206.3 


1219.6 


7 


-702.1 


-203.4 


-4305.9 


-18749.1 


8 


419.2 


125.4 


2314.3 


10789.0 


9 


-113.3 


-25.0 


-521.0 


-2650.0 


10 


68.2 


16.0 


378.0 


1715.0 


11 


-58.4 


-16.6 


-340.7 


-1533.5 



(10) 



1.2.2. Pion wave function For the 7r(140)-meson wave function ip\ ' > the solution 
obtained by spectral integral equation is rather satisfactory, it is given by the coefficients 
Ci(S — 0, L — 0, J — 0; n) which can be found in [6]. 

Still, the pion can be more precisely described by the wave function found 
phenomenologically, using the pion form factor data [37]. The phenomenological wave 
function and its parameters are as follows: 

ipn(s) = c^f exp(-&i ff s) + (3 exp(-b 27T s)^j , 

c w = 209.36GeV~ 2 , b ln = 3.57GeV~ 2 , b 2v = 0.4GeV -2 , (3 = 0.01381. (11) 
It should be noted that the difference between the wave function of Eq. (llljl and that 
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found in [6] is observed either at rather small relative momenta (k 2 = (s/4 — m?) < 0.1 
GeV 2 ) and or at very large ones. 

1.2.3. Pion emission constant The p ion-quark coupling g n for the pion emission 
q — > 7r + q (diagrams of Fig. [2] type) is given by the quark form factor g w (i4Q)-> qq -(s) 
at s = M 2 , namely, g w = g n {uo)^qq( s — M 2 ). However, the spectral integral equation 
does not determine the vertices at s < Am 2 , so in our present fit g„ is a free parameter. 

Describing the widths of p(770) — > 77r(140) and u;(780) — > 77r(140) with the use of 
vector meson ( fTUI) and pion wave functions ( fTTI) . we have found two values for the pion 
bremsstrahlung coupling: 

Solution I : \g n \ = 16.7 ± 0.3 t° 2 ;l , (12) 

Solution II :\g n \ = 3.0 ± 0.3 t° 2 i . 

The pion emission coupling, as is well known, was a subject of investigation in physics 
of low-energy pion-nucleon interactions and as well as in nuclear physics. For the pion- 

nucleon coupling, which is determined as g^NN 'n^V^I^n] , the estimations give 

g 2 NN /(<Aii) ~ 14 (see, for example, [39j EH E] and references therein). 

We can turn the description of pion-nucleon vertex into the quark language using 
quark model for nucleons: 

9nNN^P' N (Tip w ) i^lpN^j ► g^qqt^'qir^) Vytfj)^ , (13) 

see Appendix A for more detail. In Eq. ( 1T21) . we determine the vertex u — > 77T which is 
a part of the quark-language Lagrangian: 

9nqq(^P' q (r^ n ) ij B lf)^j -> 9irqq d *75^u) = 9-K d ^15^ A 

y^9nqq = #7T ■ (14) 

In Appendix A, we show that, making use of the SU(6)-symmetry for nucleons, one has 
9wNN = {5/3)9nqq- So, the SU(6)-symmetry provides us with g^-NN = (5/S^/2)g w . It 
means that Solution I does not contradict the value g 2 NN /{A'n) ~ 14 [39l HOI HI] , thus 
giving us 

16.4 < gl NN /(4n) < 23.2 . 

Note that in ( fl2l) we have included systematical errors which are due to uncertainties 
in the reconstruction of wave functions in the fit [6j. 

2. Gamma— pion decays of vector states V — > '-fir 

Here we present formulae which are used below for p — > jtt and uj — ► 77r decays. 
2.1. Polarisation vectors, amplitude and partial width for decays 

V — > 77T 



Let us introduce notations for the momenta and polarisation vectors and define the 
amplitudes and decay partial widths. 
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2.1.1. Polarisation vectors of the massive vector particle V and photon Polarisations 
of the vector meson, eVp, and of virtual photon, e^*\ are the transverse vectors: 

6 { Ppp = , e%\a = , (15) 

where q is the virtual photon four-momentum (q 2 ^ 0) and p is that of the vector 
meson (p 2 = M v ). Polarisation of the vector meson obeys the completeness condition 
as follows: 

\ " AV)( n \jy)+( n \ _ n ±v — ±p ±p _ _ iW / 1fi \ 
~ 1^ e n W e n' V a ) ~ 5V =9^' > 9p.pi ~ 5V ' ( ib J 

0=1,2,3 P 



where g is the metric tensor operating in the space orthogonal to the momentum p. 

For virtual photon, (q 2 ^ 0), the completeness condition for polarisation vectors is 
written in three-dimensional space: 

- e a*\ a ) ei S )+ ( a ) = 9a~2> , 9t2* = 9al> = 9aa> ~ • (17) 
0=1,2,3 1 

The polarisation vector of the real photon (q 2 = 0) denoted as e 7 has two independent 
components only, they are orthogonal to the reaction plane: 

= , e^p a = . (18) 
Likewise, the completeness condition for the real photon reads: 

- E # } (<0# )+ («) = g£ , (19) 

0=1,2 

„_L_L _ „ PaPa' QaQa' ± LV „ 

(pq) 

9 act' ~ 9aa' 9 n ' la — 9 oca/la' ~ la 9 Pa ■ 

2.1.2. Amplitude for the decay V —* •yn The decay amplitude V — > 777 is written as a 
product of the spin structure and form factor: 

Ag-^) = e Sg-^Cp, g)F v ^(0, M 2 ) , (20) 

with 

^ap ^ (PiO) = ^appq = ^otlWiVsPviQi^ • (21) 

In (I20j) . the electron charge e is singled out, and in (I2T]) the tensor e aia , lV2 is the wholly 
antisymmetrical. Let us emphasise the specific role of the spin operator e a p Pq . Since 
£appp — 0, this spin operator is valid for the reaction with both real (7) and virtual 
(7*) photons, so Eq. (120]) can be used for the transition with virtual photon, with 
corresponding substitution: F V ^ J7T (0) — > F^ 77r (g 2 ). 

2.1.3. Partial width for V — > 77? The partial width for the decay V — > 777 is determined 
as follows: 

2 



M v T v ^ /n = - J d§ 2 {p; q,p*) 



ap 

ap 
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^(oX)| 2 , 



a (My - Mlf 
24 Ml 



d*2(p;q,pj = ~ (9 %- jSiL (P ~ g ~ P*) • (22) 

2 (27r) d 2g (27r) d 2p 7r0 

The summation is carried out over the photon and vector meson polarisation s, and 
(e am ) 2 = (My - M 2 ) 2 /2. In the final expression a = e 2 /4ir = 1/137. 

2.2. Double spectral integral representation of the triangle diagrams with photon 
emission 

To derive double spectral integral for the form factors with photon emission by quark 
and antiquark, F^ L ^ 17T (q 2 ) and F^ L ^" I7T (q 2 ) , see Fig. [TJ one needs to calculate the 
double discontinuities of the triangle diagrams. 

2.2.1. Double discontinuities of the triangle diagrams First, consider the photon 
emission by quark, see Fig. [T^i. Corresponding cuttings for the calculation of double 
discontinuity are shown in Fig. [lb. 

In the dispersion representation, the invariant energy in the intermediate state 
differs from that in the initial and final states. Because of that, at the double 
discontinuity P ^ p and P' ^ p n . The following requirements are imposed on the 
momenta shown in the diagram of Fig. [Tb [2"3"| [55] : 

(k x + k 2 ) 2 = P 2 = s>4m 2 , (k[ + k 2 ) 2 = P' 2 =s'> Am 2 . (23) 

The momentum squared of the photon, q 2 , is fixed: 

(p - prf = (P- P') 2 = (h - k[) 2 = q 2 . (24) 

When cutting Feynman diagram, the propagators should be substituted by the residues 
in the poles. This is equivalent to the replacement as follows: (m 2 — k\)~ l — > 5(m 2 — kf), 
(m 2 — k^)^ 1 — > 5(m 2 — k 2 ) and (m 2 — k' 2 )" 1 — > 5(m 2 — k' 2 ), so the intermediate-state 
quarks are mass-on-shell: 

k\ = k\ = k' 2 = m 2 . (25) 

Then, for the diagram with photon emitted by quark (Fig. [TK), the double discontinuity 
of the amplitude (Fig. [lb) becomes proportional to the three factors: 

disc s disc^ L ^(A^) ~ Zl^G v{L) (s)G«(s') 
x d$ 2 (P; k u k 2 )d$ 2 (P'; k[, k' 2 )(27if2k 20 6 3 (k' 2 - k 2 ) 
x Sp [Q^ (L) (k) (k + m)Q^ (k[ + m)QW (-k 2 + m)] . (26) 

The first factor in the right-hand side of f[26l) consists of the following vertices: the quark 
charge factor Z^ 7n as well as transition vertices V(L) — > qq and n — > qq which are 
denoted as GV(£)(s) and G n (s'). 

The second factor contains space volumes of the two-particle states, d§ 2 (P; k%, k 2 ) 
and d§ 2 (P'; k[, k' 2 ), which correspond to two cuts shown in the diagram of Fig. [lb (the 
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space volume is determined in (I22p ). The factor (27r) 3 2k 20 5 3 (k 2 — k 2 ) takes into account 
the fact that one quark line is cut twice. 

The third factor in ( 1261) is the trace coming from the summation over the quark 
spin states. Since the spin factor in the transition V — > qq may be of two types (with 
dominant S- or dominant -D-wave), we have the following operators for virtual photon, 
Q V ^ L \ see Eq. ©: 



,LP 



LP 

9w> In' 



Q v ^\k) 



k^k 



2„,-LP 



t 1 ft 



and for the pion: 

Q« = ^75 ■ 

Here, k = {ki — k 2 )/2 is the relative momentum of the incoming quarks, k 

\ p U — LLP _ _ULP 

For real photon, we replace: 

g£(P, P')P a , = 0, g£{P, P')P' a , = , 



P 



(27) 

(28) 
ki + k 2 , 



(29) 



where (P — P') = 0. The metric tensor g^(P,P') works in the space orthogonal to 
the intermediate state momenta: g^{P, P') = g aa > - P a P a >/P 2 - P'^ p P' a + P / P$ p P'J; P . 

Actually, for the real photon we can use simpler oprator, say, = 7^, because 
in the considered decay we should have the same result for both choices, or Q^, 
due to the spin operator structure fl2Tl) . However, here we use f l29l) to emphasise an 
important point for this type of reactions: the amplitude for transversely polarized 
photons is determined by the spectral integral with transversely polarized photons in 
the intermediate states as well. 

For the photon emission, there are two diagrams: the second one is similar to that 
of Fig. [T^l but with the emission of photon by antiquark, it is shown in Fig. [lb. The 
double discontinuity of the corresponding amplitude is determined by cuttings shown 
in Fig. EH: 

disc s disc s ^^ L )^(v 7 ) ~ Zv^(v 7 )GV(L)(s)GV(s') 



x rf$ 2 (P; fcj, k 2 )d<S> 2 (P'; k[, k' 2 )(27r) 3 2k 10 5 3 (k[ - h) 
x Sp \QY (L) (k)(k + m)Q^\-k' 2 + m)Q%X-k 2 + m) 



(30) 



Likewise, there are two traces for two transitions with photon emission by quark and 
antiquark: 

SpVjf)^*^) = -Sp [Ql iL) (k)(ki + m)Q^(k[ + m)Q^\-k 2 + m) 
Spl^iV,) = -Sp <tf w (*)(*i + m)QW(4; + m)Q%\-h + m) 



(31) 



The p — ► 77r and — > 771 decays in quark-model approach 13 

To calculate the invariant form factors ^ (0) and F^ L ^ 77r (0), we should extract 

from (13"TT) the intermediate-state spin operator: 

ft 9) = , q = P-P' . (32) 

Therefore, we have: 

Spl^iV,) = Sgr^XP, q)Sl^(s, s', q 2 ) , (33) 

where 

— O^-y [b, b ,q ) , 



(^^(P,g)) 2 
(Sp^ L ^(v,)S^(P,q) 



sW-^Mf) . (34) 



2 



(S^\P,q)] 

Taking into account the expression Sp[757 ai 7 a2 7 Q3 7 a4 ] = 4de ai a 2a . iai we obtain: 

Sl?=°^(s, s', q 2 ) = Sl^(s, s>, q 2 ) = -4m , (35) 



v 

Sl{ L ^(s,s\q*)=S V v ^(s,s>,q 

with 



m 
V2 



2 , Qss'q 2 



2m z +s 



X(s,s',q 2 ) 



X(s, s', q 2 ) = {s- s'f - 2q 2 (s + s') + q\ (36) 

The photon emission amplitude, being determined by two diagrams of Fig. [T^l and Fig. 
Hb, reads 



Zl7^ L ^(q 2 , M 2 )+Z^F^ L ^(q 2 , M 2 ) 



(37) 

while the double discontinuities of the form factors in (1371) are equal to: 

disc s disc s ,F^ H7 V,M 2 ) = G v(L) (s)G v (s') 

x d$ 2 (P; h, k 2 )d$ 2 (P'; k[, k' 2 )(2nf2k 20 6 3 (k' 2 - k 2 )S V ^ L) ^(s, s', q 2 ), 
disc,disc^F^ iH7w (g 2 ,M 2 ) = G v{L) (s)G w (s') 

x d$ 2 (P; h, k 2 )d$ 2 (P'; k[, k' 2 )(2ir) 3 2k w 5 3 (k[ - k^S^^ \s, s', g 2 ).(38) 

2.2.2. The double spectral integral for the form factors with photon emission by quark 
and antiquark The equation (1381) defines the form factor through the dispersion integral 
as follows: 

°? ds 7 rf s 'disc s disc s ,PXf H7 "(g 2 ,M 2 ) 



F^^(q 2 ,M 2 )= J±j 



7T J 7T (S-M 2 ){ S '-Ml) 

4m 2 4m 2 V V(L)J\ nJ 



F v(L)^n ( 2 m2 s 7 ds 7 ds' dis Cs disc s ,F^ L) ^(q 2 ,M 2 ) 



v 7 

4m 2 4m 2 
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We have 

F^^tf, Ml) = F^ L ^(q 2 , Ml) (40) 

at equal masses of the quark and antiquark - just this case is considered here. In ( 1391) . 
we omit subtraction terms, assuming that the convergence of (1391) is guaranteed by 
the vertices GV(l)(s) and G n (s'). Furthermore, we define the wave functions of the qq 
systems: *p V (L){s) = G v{L) (s)/(s - M 2 {L) ) and ^{s') = G m (s')/{s' - Ml). 

After integrating over the momenta in accordance with fl38|) . one can represent fl39l) 
in the following form: 

dsds' 



F V ^{q\Ml)=Fl^{q\Ml)= J ^Ml)(s)Ms') 

5l^( S ,/,g 2 ), (II) 



4m 2 



9 (ss'q 2 - m 2 X(s, s', q 2 )) c v(l)^t. , 2 , 



where Q(X) is the step-function: Q(X) = 1 at X > and Q(X) = at X < 0. 

2.2.3. Z-f actors for photon emission For the p + meson, the photon emission is 
determined by two diagrams, see Figs. UK and HJo, which give us the following charge 
factors: 

Z P a^ + =*u = \ = e d = -~ . (42) 




Figure 4. Diagrams for the determination of Z-factors in the reaction p + — ► r )ir + with 
photon emission . 

For neutral vector mesons (p°, u), we have four diagrams, see Fig. which result 
in the charge factors as follows: 

Z A7 - Z V 7 - -^{eu + e d ) - - , 

Zl^° = Z^° = \{e u -e d ) = \. (43) 

In (1431) . we use the standard flavour wave functions for (I = 1, 13 = 0) and (/ = 0, 13 = 0) 

states: p° = n° = (uu — dd) / a/2 and to = (uu + dd) / 
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(d) 




Figure 5. Diagrams for the determination of Z-factors in the reactions p° — > 77r° and 
a; — > 77T° with photon emission. 



2.2-4- Decay form factors at Q 2 = —q 2 — > To calculate the integral at small Q 2 , we 
substitute: 



s = £ + ^zQ, s' = £ - ^zQ, g 2 = -Q 2 • 
In the region Q 2 <C 4m 2 , the form factors (14 ip can be written as 



(44) 



F^^(-Q 2 ,M 2 )=F^ L ^(-Q 2 ,M 2 ) = J ^V W (S)^(E) 



x 



Am 2 



-6 



7T 



1QJA(E,z,Q 2 ) 



4), A(E,z,Q 2 ) = (z 2 + 4S)Q 2 . 



(45) 



After integrating over z and substituting £ — > s, the form factors for L = 0,2 read: 



At 



(0, M 2 ) 



V7 



(0,M 2 ) = -4m J _^(a)W( )(s) 

4m 2 



x In 



s + \/s(s — 4m 2 ) 



s — J s(s — Am 2 ) 



At 



(0, M 2 ) 



v 7 



ds 



(0,M 2 ) = -m/ v / 2 / — ^(a)W(2)(a) 



4m 2 



(2m 2 + s) In 



y/s + — 4m 2 



+ 3 v / s(s-4m 2 ) 



(46) 



y/s — V s — 4m 2 

Remind that wave functions ^V(o)(s) = V'n (^ 2 )> VV(2)(s) = ij;^-' 2 ' 1 ' (k 2 ) and Vv( s ) are 
presented in Section I. 



2.2.5. Normalisation conditions for the wave functions Vv( s ) and ipv(L=o,2)( s ) It is 
convenient to write the normalisation conditions for ip n (s) and ipv(L){s) using the charge 
form factor of a meson: 



I'cliargeiy) 1 



(47) 
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The amplitude of the charge factor is defined by the photon-emission triangle diagram 
with (qq)i n = (qq) ou t- For the pion, the amplitude takes the form: 

Ax{q) = e(p + Pn) a F charge (q 2 ) , (48) 

while F charge (q 2 ) can be calculated in the same way as the transition form factors 
considered above. The normalisation condition for pion reads: 



ds , n „ , is — Am 2 



1= J T^« s » 2s V^- (49) 

4m 2 

For vector meson V(L), the normalisation condition may be determined by averaging 
over spins of the massive vector particle, see [2j [3j H21 03] for detail. Then, the 
normalisation condition reads: 



Am 2 



OO 



4m 2 

Recall that here VV(o)(s) = ^'°'^(k 2 ) and VV(2)(s) = ^ 2 ' l \k 2 ) with A; 2 = s/A - m 2 . 

2.2.6. Vector mesons: normalisation condition in case of two-component wave functions 
In the solution found in [6], the wave functions VV(o) an d VV(2) are orthogonal to each 
other with a good accuracy. Generally, vector states may mix. Then the vector mesons 
have two-component wave functions, see (jTj), and normalisation condition reads: 



4m 2 
oo 

+ ^ / 167 2 (^0 ) ^ n) +^ ( ) ^2 ) )^ 1A1 H^)^ 1 ' 2 ' 1) (^ 
4m 2 



a- (s — 4m 2 ) 2 / s — 4m 2 
x V2 - 



6 



oo 
4m 2 



(8m 2 + s)(s — Am 2 ) 2 Is — Am 2 . 
x - s -- . (51) 

3. Double spectral integral representation of the triangle diagrams with 
pion emission 

Here, we calculate the double spectral integral for the transition form factors with 
the emission of pion by quark, 7n {0, M 2 ) (diagram of Fig. [2^) and antiquark, 

F V(L)-*y*^ ( diagram of p ig _ |2J,). 
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3. 0. 7. Double discontinuities of the triangle diagrams For the diagram of Fig. [2h, the 
cuttings are shown in Fig. (2jo, with the following notations: 

fcj-^ — — k-^ — TTt ; 

(jfcj + A; 2 ) 2 = P 2 = s > 4m 2 , (JbJ + A; 2 ) 2 = P' 2 = s' > 4m 2 , 

(P - P) 2 = (*i - A*) 2 = P 2 = M 2 . (52) 

For the diagram of Fig. Eh, the double discontinuity, determined by Fig. [2b, contains 
three factors: 

disc s discy^-^(A w ) ~ Z ^/ _ >7 .(A 7r )( 77r Gy(L)(s)G 7 ( S , ) 
x d$ 2 (P; k x , k 2 )d<S> 2 (P'; k[, k' 2 ){2ir) 3 2k 20 8 3 {k' 2 - fc 2 ) 

x Sp [Q^W (fe) & + m)QW + m)g^) (-fc 2 + m)] . (53) 

The right-hand side of (133"]) is determined by the the quark charge factor Zy_ >77r (A 7r ), the 
transition vertices V(L) — > gg and 7 — > gg and p ion-quark coupling gr„.. The trace in (!53l) 
contains the operators Q^) and Q^- 1 -^ which are determined in (|29il : Qy x ^ = 7q ± (P, P') 
and Q (7r) = ry 5 . 

The diagram with the emission of pion by antiquark is shown in Fig. [2b. The 
double discontinuity of the corresponding amplitude, Fig. [2H, is written similarly to 
([53D. We have: 

disc s disc s M2 L) ^ (Vtt) ~ ^V^77r(V7r)fi , 7rG'y(L)(-S)G' 7 (s / ) 

x rf$ 2 (P; fci, A; 2 )rf$ 2 (P / ; ^) (2tt) 3 2A: 1 o(J 3 (fc; - h) 

x Sp [Q V M (k) (h + m)Q^(-k' 2 + m)Q^\-k 2 + m)] . (54) 

Correspondingly, we have two traces for two transitions with pion emission by the quark 
and antiquark: 

^'^(A 1 ) = -Sp [Ql {L \k)(ki + m)Q^(k[ + m)Q^\-k 2 + m) 

= S { ^ n) (P, P - P f )Sli L) ^(s, s', (P - P') 2 ) , 
Spl^iV.) = "Sp [Ql (L \k){ki + m)Q^\-k' 2 + m)QW(4 2 + m) 

= s$r»\p, p - p')sif^(s, s ', (p - p') 2 ) , 

S <y^){p : P-P')= e aixP{P ^ pl) = -e^ Ppi . (55) 



Here. 



(5i^(P,P-P)) 2 

(^)^( V7r )^g;^(p,p-p'; 



(s,s',(P-Pr). (56) 



2 - w v 



As a result, we obtain: 

Sl^(s, s>, (P - P') 2 ) = ^ 0) ^( S , (P - P') 2 ) = 4m 
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S v A ^(s, s', (P - P'f) = SZ^( S , s', (P - P'f 



9 6ss'(P - P') 2 

2m 2 + s + 



m 
V2 



(57) 



X(s,s',(P-P') 2 ) 

Let us note that spin factors ^l° H7 "(s, s', (P-P') 2 ) and S^f^is, s', {P-P') 2 ) differ 
by the sign only from those for photon emission 5 , ^°^ 77r (s, s', q 2 ) and S^^ jn (s, s', q 2 ), 
given by (13"5"1) . The pion emission amplitude, considered as a function of q 2 and p 2 , is 
determined by two processes (Figs. [2^, [2b): 

^iu' 1 ' "^(A"" + V,,-) = eEafipq 



X 

with 



^P^ L) ^(g 2 ,p 2 ) + 4:^P^ L) ^(g 2 ,P 2 )l , (58) 



F^^(q 2 ,p 2 n )=F^ L ^(q 2 ,pl) (59) 
due to the equality flBTI) 

disc.discyi^ H7 »*) = disc s disc s ,P^ LH7 >',p 2 ). (60) 

5. 0. 5. The double spectral integral for the form factors with pion emission The form 
factors read: 



p V(L)-> 77T/ 2 2\ _ pV(L)-> 7 7r/ 2 2\ 



rfs f ds / disc s disc s 'P^ L) ^ 77r (s',p 2 ) 



7T 7 7T (s-M 2 fL 0(s'-g 

W 4m 2 V 



(61) 



As in ( 13"9"j) . we assume that the convergence of flBTT) is guaranteed by the vertices Gym(s) 



and G 7 (s'). 



Futhermore, we consider the production of photon, q 2 = 0, and use the photon wave 
function t/> 7 (s') = G 1 {s')/s' . After integrating over intermediate-state quark momenta, 
one can represent f lBTl) for p 2 < in the following form: 

^A" ( U )PJ - l U >PJ 



oo 



4m 2 

x Q^^-^'^^^^^^^p 2 ). (62) 

V A 0> s '^) 



The step-function 0(X) was defined in (jHJ). 

Let us emphasise once again that Eq. (I6"2"j) is valid in the region p 2 < only. 
To obtain form factors at p 2 = M 2 , one needs to continue Eq. (1621) to the region 
p 2 > 0. Since the form factors are analytical functions in the vicinity of p 2 = 0, the 
straightforward way is to expand them in a series over p 2 keeping constant and linear 
terms only: 

P^ L) ^(0,p 2 ) = Fl^(0,0) +Pl-^FZ^ (0,0) . (63) 
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One can approximate pi -^^^(O, 0)/dpl = ^^^(O, 0)-F$ iH77r (0, -Ml) (here 
^ = -Ml). Then 

i^ L) ^(0,M 2 ) = ^'^(0,^) 

= 2F^ LH7 "(0, 0) - FS L) -^(0, -M 2 ) . (64) 

Both form factors, "^(O, 0) and F^ i ' ) ^ 77r (0, — M 2 ), are calculated according to 
Eq. (1521). 



3.0.9. Z-f actors for pion emission The charge factors for the pion emission in the 
decays p + — > 7vr + , p° — > jtt°, uj — > 77?° (see Figs. El [7J are equal to those for photon 
emission as follows: 

/yp^ — >77r^ — ^7 7r ^ — >77r~*~ — ^T 71 " - *" 

7 P°-7-° _ 7 P°-7-° _ 7 P°-7-° _ 7 P°-7-° _ 1 ( p _\_ p \ 

Z17^° = = = Z£T*° = \(e u - e d ). (65) 

In the calculation of Z-factors (|65|) . we take into account that pion emission by quark 
is a two-step process (see Fig. [6b ): the initial quark (for example, in Fig. [6k.) emits 
gluons (they have isospin I g i uon s = 0) which produce quark-ant iquark pairs, uu or dd, 
with equal amplitudes, and then we face the transition ud — ► 7r + . The block of Fig. [6b 
is denoted as a coupling 

In the process of Fig. [Th,, the gluons produce uu pair with the same amplitude as in 
the previous case but then we face the transiton uu — > tt° resulting in the factor l/y/2 
(recall that tt° = (uu — dd)/y/2). In the process of Fig. [7b, the dd pair is produced, and 
the transiton dd — > n° gives the factor —l/y/2 (for more detailed presentation of the 
quark combinatorial rules see [2] and references therein). 




Figure 6. Diagrams for Z-factors in the reaction p + — » jn + : a) Z — ed — — | and 
b) Z = e u = §. 



3.0.10. Partial width In terms of the calculated form factors, the partial width reads: 

1 aM*-Ml\(M*,Ml,Q) 
3 ' 4 



M V T 



V 1 - V->~{1T 



X 



4 My 2 

jXr^r^o, m 2 ) + zXr 7w Fir^(o, m 2 ) 



+Z^Fl^(0, Ml) + Z^F^(0, Ml) 



(66) 




Figure 7. Diagrams for Z-factors in the reactions p° — > jir° and ui° — > -fir : a) 

Z(p° -> 7 ^°) = i, Z( W - 7^°) = b b)Z(p° -> 7 ^°) = |, - 7^°) = f; 

c)Z(p° -> 7 7r°) = -|, -> 7^°) = i; d)Z(p° -> 7 ^°) = -I -> 7 ^°) = i 
Recall that p° = tt° = Hii^M, and uj° = a^+M. 



Here, the factor 1/3 is due to the averaging over initial vector meson spin states, the 
term 

a/4 • (My — Ml) /My is given by the phase space integration, and A(My,M^,0)/2 = 
(My — M^) 2 /2 is due to the spin factor (|2~T|) . The Z-factors are as follows: 
z£T*° = 1/6, = 1/6, z£^° = 1/6, Z^T* = 1/6, 

= 1/2, = 1/2, = 1/2, Z^f = 1/2. 

4. Results and discussion 

The fitting to the partial widths T ( p ±^ n± = 68 ± 30 keV, T^^o = 77 ± 28 keV, 
r^^Lo = 776 ± 45 keV leads to the following values of the pion emission coupling: 
Solution I : 16.7 ±0.3 t° 2 i , 

Solution II: - 3.0 ± 0.3 t 2 \ ■ (67) 

In Eq. ( ISTl) . we have included systematical errors ((+0.1/ — 2.3) for Solution I and 
(+2.1/ — 0.1) for Solution II) which are caused by the uncertainties of the fit of qq wave 
functions in the spectral integral equation (see Section 1.2). 

So, we have regions of positive and negative g n . However, one should take into 
account that the sign of g n in (JBTl) is rather conventional: it depends on signs of wave 
functions of photon and mesons involved into calculation. Because of that, being precise, 
we should state that for g n we determine absolute values only, see ffT2l . 

Solution I gives us the value of the of pion-nucleon coupling; recall that it is 
determined as a factor in the phenomenological Lagrangian: g^NNi^'ni^^l^N^)- 
It is in agreement with the results for pion-nucleon scattering gl NN /4iT ~ 14 [3"9 | |4"0"1 |4"T] . 
Namely, dealing with pion-nucleon interaction in terms of the quark model, we use the 
Lagrangian: 

g*qq($'q(T<PTr) «75^g) = ^9irqq V>t+ ($ 'd + Other terms 

= N> ' d i^u^j + other terms, (68) 

that gives us y/2g„ w = </ w . 
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In Appendix A, using SU(6)-symmetry for nucleons, we demonstrate that g n NN = 
{h/3)g 7:qq . So, in terms of SU(6)-symmetry, we have: 

5 . . 

9nNN = ^^9n- (69) 



We see that Solution I, being in agreement with data [3H1 SOI EH] , gives us 

gl NN /(4n)= 22.2 ± 0.8±°;g. (70) 

For Solution II, we have found 0.03 < g^ NN /{A'K) < 1, that is far from the 
experimental value. 

4-1. Predictions for excited vector states 

For p ± (2S) , p°(2S) and u(2S) mesons, we have found the following partial widths (in 
keV units): 

r(/4^ 7 0-io-i3o, 

r(p^^7vr)-10-130, 

r(cj 2 5-*77r) - 60 - 1080. (71) 

The other wave functions of highly exited states have too large uncertainties to provide 
us with reliable widths. This points to the necessity to carry out mesurements of 
radiative processes with mesons in the region of large masses. 

Acknowledgement We thank B.L. Birbrair for helpful remarks. This paper was 
supported by the RFFI grant 07-02-01 196-a. 

Appendix A: Nucleon pion emission vertex 
in the SU(6) quark model 

Here we derive the relations between couplings in phenomenological Lagrangian for 

pions and nucleons, gnNNypNi^'f^l^N^j, and those for quarks, g vqq (^' q {T(p 7 ,)i^ 5 'ip^j. 

To be definite, we consider transitions p^ — > n + + and — > ir + + dK We use the 
following SU(6) wave functions (see, for example, Appendix D in Ref. [44]): 

ip p = p T (q(l)q(2)q(3)) = ^(wVd 1 + dVw T + w T dV) 

{v)u^d) + d)v)u^ + u^dJv) + u^u^d) + d)u^v) + v)d)u^). 



3^/2 

V2 



= n l (q(l)q(2)q(3)\ = ^-(d l d l u 1 + u^d l d l + d l uU l ) (72) 
L -.(d l d^u l + u l d l d^ + d^u l d l + d^d l u l + u l dU l + dVd T ). 



3y/2 

Recall that for baryon quarks we use notation of the type d^u^d) = £^(1)^(2)6^(3). 
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The isospin block reads: 

Transition — > 7r + + is given by the following terms in nucleon and quark spaces: 

ftrAwfa^l V2 175b 1 ) 

= ^ M <7T + n^g(l)g(2)g(3))| V2 

xE 3 ^y^ ^(ffdM2)ff(3))> , (74) 

where n(j), T 2 (j) and 75 (j) act on 

In the non-relativistic limit, which we use for nucleons and constituent quarks, 



*7s -> (-i)(Zq) 



and direct calculations give: 

5 

flV/ViV = S^gg * 3 • - . (75) 

To simplify the calculations which lead to (175]) . one can fix the direction of photon 
momentum, for example, q = (q x ,0,0) and then use (qV{j)\ — > qx&iij)- 
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